Path derivation for a wave scattered model to estimate height correlation function of 

rough surfaces 
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The long standing problem on finding the height correlation function is studied by the inverse 
scattering problem. We propose a new method in the frame work of Kirchhoff theory which we call 
"path derivation of scattered wave (PDSW)" in order to obtain an expression for direct measure- 
ments of the height correlation function. This would provide adequate insight to rough surfaces. 
The efficiency of this method is due to the fact that the height correlation function could be esti- 
mated directly by measurements of the scattered intensity on a suggested path. The model is tested 
numerically and an experimental setup is suggested. 

PACS: 42.25.Fx, 02.30.Zz, 02.50.-r 



CZ5 

o 

H— > . 

Or 

q 
o 

Or 



> 

O 
O 

cn 

o 



x 



PACS numbers: 

The importance of rough surfaces in various fields of 
science and engineering has been highlighted by many 
researchers. The area of study spans from Biology [l|, 
Polymer Q and Radio Science, Material engineering 
[3], and Physics (see Q and references therein), etc. This 
motivates researchers specially in Physics to study the 
various aspects rough surfaces. The stick slip dynamics 
of rough surfaces was studied experimentally for two 
surfaces where a hydrodynamic model was implemented 
to theoretically study the effects of surface roughness on 
viscous incompressible liquids Q- Palasantzas studied 
the effect of roughness on ballistic thermal conductance 
of a nanosized beam in addition to effects of the rough- 
ness exponent on friction while a viscous rubber slips in 
to a rough surface In a detailed study Volokitin and 
Persson [lfj| stated that if surface modes for instants ad- 
sorbate vibration modes exist on a surface, thermal radi- 
ation would be coherent and due to the electromagnetic 
oscillations on the surface, radiative heat transfer would 
be experienced. 

The topography of any surface could be obtained by 
mechanical devices such as scanning probe microscopy, 
SPM [ll|, [l2| e.g. in application to semiconducting car- 



bon nanotube transistors [131 ] . But despite the wide range 
of applicability and benefits of SPM, the tip convolution 
may still remain as the main challenge for SPM (lij . in 
addition to the limitations that SPM has in changing 
the scale of scanning. Optical techniques would over- 
come this limitation without scratching the surface. Op- 
tical techniques could also manage to determine the scale 
of observation with the incident light wave-length. The 
mostly implemented method in the optical technique is 
the scattering method where the roughness would play 
a major role on the scattering of waves [Hf. Interest- 
ing work in the context of wave scattering has been car- 
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ried out in application to Brewster's scattering angle [16( 
and the Rayleigh hypothesis [13] ■ Since the wave-length 
could be either smaller or greater than the height fluc- 
tuations, the two scale theory was proposed [18j. In a 
further study, effects of interference of two beam scat- 
tering was considered [l||. All these plus the fact that 
studying the scattered wave would provide very useful 
information about the surface itself [20j, would motivate 
studying the inverse scattered problem. 

Inverse scattering techniques were implemented 
[iij l to measure the statistical properties of rough surfaces. 
Although after quite a long time since optical techniques 
were born, and many parameters like; root mean square 
of height [2^|, probability density function (PDF) [27j |. 
correlation length [28| etc. where estimated, not much 
progress has been achieved in order to find the most im- 
portant parameter of the surface which is the height cor- 
relation function. This issue is the main aim of this work. 
The framework implemented is wave scattering by Kirch- 
hoff theory approach [TBI, l29l - [3l| . The Kirchhoff theory 
is an electromagnetic theory which treats any point on a 
scattering surface as a part of an infinite plane, parallel 
to the local surface tangent [l5| . In this work a model 
is proposed to estimate the height correlation function 
of a random surface which we name; " path derivation of 
scattered wave" . Taking into account this model enables 
experimentalists to obtain the height correlation function 
simply by only measuring the intensity of the scattered 
light in a special path. This technique would dramat- 
ically simplify the estimation of the height correlation 
function due to the fact that only one parameter (the 
scattered intensity) needs to be measured. 

In Kirchhoff theory, the field of incident monochro- 
matic wave may be written as tp mc (r) = exp(— ifc.;„ c .r), 
where k and r indicate the wave number and position re- 



2 



Incident 
beam 




Scattered 
beam 



FIG. 1: The geometry used for wave scattering from a rough 
surface. 



spectively. In Fig. 1 a beam is shown incident to a rough 
surface with dimensions — X < xo < X , — Y < yo < Y 
with a reflection coefficient Rq equal to —1 corresponding 
to the Dirichlet boundary condition. The total scattered 
field over the mean reference plane Am, is given by [l5| 
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sm t>i — sin 6>2 cos 0, 

— sin 02 sin 0, 

— (cos 6>i + cos 62), 

sin 0i (1 — Rq) + sin 62 cos 0(1 + Ro), 
sin 6*2 sin 0(1 + Rq), 

cos6» 2 (l + .R )-cos(9i(l-i? ). (2) 



The total scattered intensity for isotropic surfaces could 
be considered as sum of the diffused and coherent parts, 
the diffused part is [15] 
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where C'or(R) = W X +KW X )) j s t ne height- height corre- 
lation function and (h(x)) — 0, g — fc 2 <7 2 (cos#i + cos#2) 2 
is the roughness criteria, and a — -J ((/i(x) — (ft,(x))) 2 ) 
is height root means squared of the rough surface. It is 
more convenient to write equation ([3]) in the form 
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with Id/ Id 



^£rA M exp(-g), 



where from hereafter the 



overtilde is omitted. Equation (j4|) shows the diffused in- 
tensity in terms of the height correlation function, where 
the idea is to find the height correlation function in terms 
of the diffused intensity. This would be obtained by solv- 
ing equation ((4]) making use of the saddle point approx- 
imation. Attention must be focused on a few issues; we 
are considering rough surfaces were 5 > 1, this means 
that the conditions, ka ^> 1 and small 0\ and 62 must be 
fulfilled. 

- In highly rough surfaces, g > 1, the contribution 
of the coherent term could be neglected in comparison to 
the contribution of the diffused term, this means that the 
diffuse intensity Id would be equal to the total intensity 
I [HI, expressed as 

I = J Jo (kRy/A 2 + B 2 j exp(gCor(R))RdR. (5) 

- In order to satisfy the conditions of the saddle point 
approximation, the integrand must tend to zero on the 
integral limits. This condition is not satisfied by Eq. (|5|), 
persuading us to differente both side of Eq. ([5]) in respect 
to g resulting in 
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Note that the derivative has been taken on a path where 
the Bessel function Jo(fci?V A 2 + B 2 ) stays constant in 
respect to the variable g. It could readily be seen that 
the term Jo(/ci?V A 2 + B 2 ) depends on the wave-length, 
incident and scattered angles. The Bessel function is kept 
constant with respect to the variations of g. This is done 
by keeping the wave-length and incident angle constant 
and only letting the scattered angles vary in a way to 
keep the A 2 + B 2 constant, see Fig. 2. In Fig. 2, three 
typical paths corresponding to three different constant 
values of A 2 + B 2 is shown. 

Making use of the definition of the logarithmic function 
we obtain 
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Equation would satisfy the conditions for the saddle 
point approximation, where the argument of the Bessel 
function must stay constant in respect to g. According 
to Fig. 2, the variation of 82 and would provide the 
path where the intensity is measured on. Note that in 
plotting the curves (paths) of Fig. 2, the Bessel function 
stays constant with respect to the variations of g. 

- To solve the integral of Eq. ((7|) by the saddle point 
approximation, Eq. ([7]) needs to be written in a form 
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FIG. 2: (Color online) The curves are typical paths on which 
the Bessel function or A 2 + B 2 stays constant with fixed k as 
g varies. The plane of consideration is the diffused intensity 
based on the scattered angles 82 and <j> for a Gaussian corre- 
lation function. The values for the parameters are; 9\ = 20°, 
a = 50nm and A = 500nm. The paths correspond to (a) 
A 2 + B 2 — 0.7 (b) A 2 + B 2 — 1 (c) A 2 + B 2 = 1.2. 



so how that the extremum would be inside the integral 
limits and not on the integral limits. By implementing 
a new variable u = hi(R), the lower limit of the integral 
could be changed to —00. Hence, Eq.fJT]) would be 
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Since the intensity and intensity variations are needed on 
the path, and depending on the chosen constant A 2 + B 2 , 
one could refine the correlation. This could be obtained 
due to the fact that all paths would give a unique height 
correlation function. This would suggest a new technique 
for experimentalists to test their obtained results by 
choosing specific incident angles, incident wave-lengths, 
to keep the desired Jo constant in respect to g. In other 
words to keep Jo constant, the parameter A 2 + B 2 must 
stay constant. In this way a new path may be suggested 
for finding the intensity experimentally. Fig. 2 shows 
three paths obtained for three constant values of A 2 + B 2 , 
6»i = 20° and A = 500 nm. 

In this stage Eq. |[8]) complies the saddle point approx- 
imation by 
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where, u* is the desired saddle point. Solving the integral 
in equation ([8J by the saddle point approximation gives 
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(gf"(u*))-hM9f(u*)), 



(10) 



due to the complexity of the RHS of Eq. (JT0J) , an explicit 
version for the height correlation function may not be 
obtained without the use of the second derivative of Eq. 
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Equation (|11[) is an expression in terms of Cor(u) in 
the extremum point which itself is dependent on g. This 
means that changing g by differing the wave-length, stan- 
dard deviation, incident and scattering angles, the cor- 
relation function would vary. This enables us to obtain 
the functionality of the correlation on g instead of just 
having a specific point. By dividing both sides of Eqs. 
(fTTJl) and (ITT|) . the expression for Cor(u*) is obtained 
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(12) 



Equation ([12")) shows the relation between Cor(u*) and 
g. Due to the fact that W(g) could be observed exper- 
imentally, Eq. (TT2|) would give the correlation function. 
But the issue which still remains is the value of u* . 

From the second expression of Eq. ((9]) an explicit rela- 
tion between the derivative of the height correlation and 
g for the limiting case of highly rough surfaces (j > 1) 
is obtained 



Cor (u*) = — 
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The derivative of Eq. (fT2"j) with respect to g would 
show the functionality of u* in terms of g 



dCorju*) du*~* dW{g) 
dg dg dg 



and by using Eq. (13) 
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which with partial integration, it could be obtained u* (g) 



u*(g) 



9W{g) 



W(g) 



dg. 



(16) 



Equation (1161) gives an expression for u* on g. The re- 
lation between Cor(u*) with g is already shown in Eq. 



4 




1 1.5 2 2.5 3 

R 



FIG. 3: Comparison of the correlation function estimated by 
the PDSW model and the original exponential correlations. 
The parameters for this estimation are; A 2 + B 2 — 0.8, 9 1 — 
20°, a = 0.5 and ka = 1.2. 

(fE2")l . so the dependence of the height correlation func- 
tion on u may readily be obtained. Since u = ln(i?), the 
dependence of the height correlation function on R may 
be obtained. 

The procedure could be illustrated as follows; inci- 
dent light with a specific angle to a rough surface with 
a constant wave length (preferably small) scatters with 
a verity of angles. The measured intensity for each scat- 
tered angle on the desired path needs to be differentiated 
in respect to g (which depends on the scattered angle) 
twice. Then the ratio of the second derivative and the 
first derivative is obtained in order to find W(g). Substi- 
tuting W(g) in Eq. (fT()|) the parameter u*(g) is obtained, 
keeping in mind that W(g) is the height correlation of 
u* . ^Having the dependence or in other words variations 
of Cor(u*) and u* on g the height correlation on R is 
obtained. 

In Fig. 3 the applicability and robustness of this tech- 
nique is illustrated. Where as an example we estimate the 
known exponential correlation function with the PDSW 
model and compare it with the original cases showing a 
very good consistency. The parameters for this estima- 
tion are; A 2 + B 2 = 0.8, = 20°, a = 0.5 and ka = 1.2. 

The applicability of the method introduced in this 
work for experimentalists although has a very simple pro- 
cedure, would give the height correlation function of the 
surface. The height correlation function is the most im- 
portant feature in order to understand the physics of 
a rough surface which is applicable to various fields of 
physics. "The path derivative scattered wave" intro- 
duced in this work would prove as an easy to use tech- 
nique for experimentalists since only by measuring the in- 
tensity of the scattered light on a specific path, the height 
correlation function of a specific rough surface could be 
obtained. 
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